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摘 要: 养殖 场 中 肉牛 较为 活跃 ， 采 集 得 到 的 图 像 数据 中 肉牛 姿态 多 变 ， 肉 牛 姿态 端正 帧 较 少 ， 导 致 自动 
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骨架 特征 和 肉牛 图 像 边缘 轮廓 特征 ， 提 出 一 种 多 


姿态 肉牛 体 尺 自动 测量 方法 。 首 先 ， 利 用 深度 相机 Azure Kinect DK 从 正 上 方 采集 肉牛 俯视 深度 视频 数据 ， 
对 视频 数据 进行 分 帧 处 理 ; 其次， 对 原始 深度 图 像 进行 预 处 理 ， 将 肉牛 从 复杂 的 背景 中 提取 出 来 ; 再 次 ， 


利用 Zhang-Suen 算 法 提取 目标 图 像 肉 牛 骨 架 ， 检测 


骨架 交点 和 端点 ,分 析 肉 牛头 部 特征 ， 并 确定 头 部 去 除 


点 ， 去 除 图 像 中 肉牛 头 部 信息 ; 最 后 ， 利 用 改进 的 U 弦 长 曲率 算法 提取 肉牛 轮廓 曲率 曲线 ,根据 曲率 值 确定 
体 尺 测 点 ,将 体 尺 测 点 转换 到 三 维 空间 中 ,计算 体 尺 参数 。 本 研究 通过 分 析 大 量 深 度 图 像 数 据 ， 将 图 像 中 
肉牛 姿态 分 为 左 牌 、 右 牌 、 姿 态 端正 、 低 头 和 抬头 五 类 。 试 验 结果 表明 ， 本 研究 提出 的 基于 骨架 的 多 姿态 
肉牛 头 部 去 除 方法 在 5 种 姿态 下 的 头 部 去 除 成 功率 均 高 于 92%; 在 23 头 肉牛 不 同 姿态 共 46 帧 深度 图 像 中 ， 


利用 基于 改进 U 弦 长 曲率 的 体 尺 测 点 提取 方法 ， 测 得 体 直 长 测量 的 平均 绝对 误差 为 2.73 cm， 体 高 测量 的 平 
均 绝 对 误差 为 2.07 cm， 腹 宽 测量 的 平均 绝对 误差 为 1.47 cm。 研究 结果 可 为 精确 测量 多 姿态 下 肉牛 体 斥 提供 


支撑 。 
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1 引 Ë 


随 着 中 国人 民 腾 食 结构 的 改善 ， 人 们 对 牛肉 
的 数量 和 品质 要 求 愈 来 愈 高 ， 联 合 国 粮 食 及 农业 
组 织 最 新 数据 显示 ，1990 年 至 2020 年 中 国人 均 
牛肉 的 消费 量 上 涨 5.25 倍 …。 肉 牛 体 尺 参 数 是 肉 
牛 生长 发 育 状 况 以 及 选 种 育种 的 重要 判断 依据 ， 
而 在 现 阶 段 ， 大 多 数 养殖 场 仍 采 用 传统 人 工 方法 
测量 肉牛 体 尺 数据 ”， 这 种 测量 方法 工作 量 大 、 
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效率 和 准确 性 低 、 容 易 造 成 肉牛 应 激 ， 还 会 影响 
动物 福利 ， 因 此 开展 肉牛 体 尺 的 自动 化 测量 研究 
具有 重要 意义 =。 

国内 外 已 有 较 多 基于 计算 机 视觉 的 牲畜 体 尺 
测量 技术 的 研究 ， 且 取得 了 较 好 的 研究 进展 。 
Spoliansky 等 ”使 用 Kinect 相 机 自动 提取 牧场 中 
每 头 奶牛 的 个 体 特 征 ， 并 设计 回归 模型 对 奶牛 的 
身体 状况 评分 ， 结 果 表明 机 器 学 习 方 法 可 自动 评 
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测 奶 牛 体 况 。Nir 等 :使 用 单个 Kinect V2 传感器 
获取 奶牛 深度 图 像 ， 采 用 机 器 学 习 方法 、 椭 圆 拟 
合 方法 和 分 位 数 回归 方法 提取 奶牛 肩 高 、 辟 高 等 
体 尺 参数 ， 结 果 表 明 利 用 深度 相机 采集 到 的 深度 
信息 可 自动 测量 奶牛 体 尺 参数 ， 且 测量 精度 较 
高 。Salau 等 使 用 六 个 Kinect 相 机 采集 奶牛 深 
度 图 像 ， 结 合 手动 标记 的 前 乳头 和 坐骨 结 节点 ， 
通过 标记 点 成 功 计算 出 静止 和 行进 状态 下 奶牛 坐 
骨 长 度 和 坐骨 结 节 的 高 度 。 赵 建 敏 等 ”提出 一 
种 基于 Kinect v4 传 感 器 的 肉牛 体 测量 方法 ， 采 
集 彩 色 和 深度 图 像 ， 结 合 目标 检测 和 边缘 检测 等 
方法 ,测量 牛 体 嗜 甲 高 、 体 斜 长 和 体 直 长 等 参 
数 ， 体 尺 测 量 误差 仅 为 0.76%、1.68% 和 2.14%。 
李 嘉 位 等 S EJ Kinect DK 深度 相机 、 红 外 对 射 
光栅 触发 器 和 射频 识别 搭建 三 维 点 云 采 集 系 统 ， 
实现 了 肉牛 点 云 的 三 维 重 建 与 分 析 ， 采 集成 功率 
为 91.89%， 体 尺 重 建 误差 为 0.6%。 赵 新 强 ” 通 
过 凸 包 算法 去 除 牛 涉 和 牛尾 获取 奶牛 身体 区 域 ， 
通过 角 点 分 析 方 法 获取 奶牛 体高 、 体 宽 、 体 长 数 
据 ， 体 尺 平均 误差 均 在 3.2% No WER O 将 
牛 体 姿态 分 为 抬头 和 低头 ， 在 不 同 姿态 下 设计 不 
同 的 体 尺 测量 算法 ， 实 现 了 牛 只 体高 的 测量 ， 体 
高 测量 的 平均 误差 为 1.9%， 胸围 和 腹 围 的 预 估 
误差 分 别 为 2.865% 和 3.463%。 以 上 研究 中 ， 特 
畜 体 尺 测量 多 是 基于 姿态 端正 图 像 数 据 ， 但 在 实 
际 自动 化 测量 应 用 中 ， 这 些 被 测 牲畜 多 是 处 于 运 
动 状 态 ， 采 集 图 像 中 的 牲畜 姿态 也 多 是 不 端正 
的 ， 导 致 性 畜 体 尺 测 量 效率 较 低 。 

针对 上 述 问 题 ， 为 提高 牲畜 体 尺 测量 效率 ， 
本 研究 利用 肉牛 骨架 特征 和 肉牛 轮廓 边缘 特征 ， 
开展 图 像 中 多 姿态 下 肉牛 头 部 去 除 方 法 和 体 尺 测 
点 提取 方法 研究 ， 实 现 自动 测量 多 姿态 下 肉牛 体 
尺 参 数 。 


2 材料 与 方法 


2.1 数据 获取 


试验 数据 采集 于 河北 省 保定 市 定 兴 县 某 肉牛 
养殖 场 ， 采 集 对 象 为 10 月 龄 左右 的 西门 塔 尔 肉 


牛 共 34 头 ， 采 集 时 间 为 2021 年 9 月 30 日 。 数 据 
采集 区 域 长 2.37 m， 宽 1.53 m， 采 集 装 置 如 图 1 
所 示 。 当 肉牛 进入 数据 采集 区 域 后 ， 关 闭 前 后 
门 ， 利 用 深度 相机 采集 肉牛 俯视 深度 视频 数据 ， 
采集 完毕 后 打开 前 门 ， 肉 牛 进入 休息 区 域 。 深 度 
相机 支架 架设 在 墙 体外 侧 ， 以 降低 牛 的 应 激 反 
应 。 调 整 相 机 高 度 和 角度 使 其 能 拍摄 到 完整 的 肉 
牛 图 像 ， 相 机 距离 体重 秤 所 在 平面 2.92 m， 调 整 
相机 为 NFOV2X2 装 箱 模式 ， 帧 率 为 15 fso F 
用 ffmpeg 工 具 对 采集 到 的 肉牛 俯视 深度 视频 数据 
进行 分 帧 ， 得 到 肉牛 俯视 深度 图 像 。 


Er 
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(a) 数 据 采 集 装 置 示意 图 
数据 采 


oxen 。 集 区 8 TRAN ae 


Ab fe 


(c) 肉 牛 活动 区 域 示意 图 
图 1 肉牛 数据 采集 装置 


Fig. 1 Data acquisition device of beef cattle 


2.2 图 像 预 处 理 


深度 相机 采集 得 到 的 图 像 数 据 存在 一 定 畸 
变 ， 为 检测 相机 在 数据 采集 时 的 误差 ， 并 对 
图 像 进行 畸变 矫正 , 将 相机 架设 在 距离 尺寸 为 
8 cm X15 cmX29 cm 的 长 方 体 物体 前 方 1.5 m 
处 ， 利 用 Matlab 软件 标注 实际 采集 到 物体 大 小 ， 
数据 采集 误差 为 2.49 cm， 如 图 2 (a) 所 示 。 为 
减 小 数据 采集 误差 ， 本 研究 利用 8X9 标 定 板 对 
相机 的 内 部 参数 和 畸变 系数 进行 标定 ， 标 定 后 数 
据 采 集 误 差 为 0.35 cm， 如 图 2 (b) 所 示 ， 满 足 
数据 采集 需求 。 

由 于 养殖 场 环 境 复杂 ， 采 和 集 到 的 原始 深度 图 


Xiv 合 作 期 刊 


146 
0.30 o 
nap 0.25 
DaS 0.20 
oia 0.15 
0.10 a0 (a) 深 度 差 值 图 (b) 归 一 化 图 
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(a) 原 始 数据 采集 效果 (b) 矫 正 后 数据 采集 效果 
图 2 相机 矫正 前 后 肉牛 数据 采集 效果 (ce) 分 割 后 图 像 (d) 目 标 图 像 


Fig. 2 Data collection effect of beef cattle before and after 
depth camera correction 

像 中 除 包含 肉牛 外 ， 还 有 栏杆 、 地 面 和 地 磅 秤 等 

复杂 的 背景 信息 ， 这 些 背 景 信息 影响 体 斥 测 点 的 

提取 精度 ， 因 此 本 研究 去 除了 图 像 中 的 背景 信 

息 ， 处 理 效 果 如 图 3 所 示 ， 具 体 步 骤 如 下 。 

(1) 利用 影 差 法 '" 获取 深度 图 像 f(x, y) 
与 背景 深度 图 像 f(x, y) 的 深度 差 值 图 像 
D(x, y)， 计 算 如 公式 (1) 所 示 ， 
图 3 (a)。 

D(x,y) = f(x,y) -f(x,y) (1) 

本 研究 中 青绿 色 图 像 均 是 经 处 理 后 的 可 视 化 
图 像 ， 原 始 图 像 为 16 位 深度 图 像 ， 可 视 化 效果 
不 佳 。 

(2) 对 16 位 深度 图 像 进 行 归 一 化 "” 处 理 ， 
得 到 图 3 (b) 所 示 图 像 。 图 3 (a) 中 像素 最 高 
EN D na Yn) RR REEN D Onn Ynn) 
则 归 一 化 后 (x, y) 点 的 像素 值 D(x, y) 变 为 
Drom WAR (2)。 

D (x,y) -D (Xunn Ynn) 
D (Xna Ymar) — D (X mins Yin ) 

(3) 利用 基于 分 水 岭 的 图 像 分 制 ”方法 获 
取 含 有 牛 的 二 值 图 像 ， 但 由 于 牛 与 栏杆 相距 较 近 
且 深 上 度 值 差 异 不 大 ， 导 致 分 割 后 的 图 像 带 有 少量 
栏杆 像素 点 ， 如 图 3 (c) 所 示 。 

(4) 利用 提取 图 像 连 通 域 方 法 得 到 只 含有 肉 
牛 的 目标 图 像 ， 如 图 3 (d) Bras. 


结果 见 


用 页 = x 255 (2) 


图 3 肉牛 深度 图 像 预 处 理 结果 


Fig.3 Deep image preprocessing results of beef cattle 
2.3 基于 骨架 的 多 姿态 下 肉牛 头 部 去 除 


针对 牛头 摆动 造成 体 尺 测量 误差 较 大 的 问 
题 ， 赵 新 强 采用 凸 包 分 析 方 法 去 除 牛 头 信息 ， 
将 牛 只 身体 区 域 作为 体 尺 测 量 对 象 。 但 凸 包 分 析 
方法 只 适用 于 牛 姿态 端正 图 像 ， 鲁 棱 性 较 差 。 

为 去 除 多 姿态 下 的 肉牛 头 部 信息 ， 本 人 研究 提 
取 了 肉牛 背部 主 骨架 ,效果 如 图 4 所 示 。 首 先 ， 
利用 Zhang-Suen 算 法 ”不断 和 迭代 细 化 处 理 肉 牛 
目标 图 像 C(x，y)， 得 到 肉牛 原始 骨架 ; 其 次 ， 
利用 单 像素 细 化 算法 将 肉牛 骨架 宽度 细 化 为 1 像 
A; 再 次 ,利用 3X3 卷 积 核 遍 历 图 像 中 的 前 景 
点 ， 提 取 骨 架 交 点 和 端点 ， 并 设置 骨架 交点 处 8 
邻 域内 像素 值 为 0， 通 过 最 大 连通 域 方法 将 肉牛 
整体 骨架 分 为 若干 个 骨架 ， 如 图 4 (a) 所 示 ; 最 
后 ， 计 算 图 4 (a) 中 各 骨架 最 小 外 接 和 矩形 所 占 面 
只 ,得 到 和 矩形 面积 最 大 的 骨架 即 为 肉牛 主 骨架 ， 
如 图 4 (b) 所 示 。 

交点 约束 条 件 为 公式 (3)。 


os a (3) 
S(P) = 255 

端点 约束 条 件 为 公式 (4)。 

TEA (4) 
S(P)= 255 


其 中 ,N(P) 是 与 P 相 邻 的 8 个 像素 中 的 前 景 
点 个 数 ， 个 ; S(P) 为 P 点 处 的 像素 值 。 
扫描 图 4 (b) 中 主 骨架 像素 点 ， 存 在 两 个 8 
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(a) 骨 架 交 点 与 端点 位 置 (b) 主 骨架 
图 4 利用 Zhang-Suen 算 法 获取 肉牛 主 骨 架 
Fig. 4 Getting the main skeleton of beef cattle by 


using Zhang-Suen algorithm 


邻 域内 N(P)= 1 的 点 。 比 较 两 点 间 坐 标 位 置 关 
系 ， 位 置 偏 左 靠 上 的 点 为 主 骨 架 前 端点 P,， 男 一 
点 为 主 骨 架 后 端点 P,。 

由 于 图 像 中 肉牛 骨架 为 离散 点 ， 利 用 后 向 差 
分 方法 求 忆 ,点 处 的 斜率 。 设 骨架 前 端点 为 
G, (xp yi) A G1 为 起 始点 依次 计算 主 骨架 序列 
G(x, yi) 中 相 邻 两 点 间 的 距离 4， 用 NN 记录 4d = 
V2 出 现 的 次 数 , 若 X= 3 则 记录 最 后 一 次 参与 计 
SUPRA G(x, yi) ARR. EAR BiG 


点 处 的 斜率 计算 如 公式 (5)。 
o (5) 
分 割 线 的 直线 方程 为 公式 (6)。 
yao (xox) ty, (6) 


利用 Canny 算 子 提取 C(x, y) ARK fe 
B, AHHAR (6) 计算 得 到 直线 ， H 
分 为 两 部 分 ， 比 较 两 部 分 面 只 ,保留 面积 较 大 区 
Ih A(x, y)。 

利用 凸 包 分 析 的 方法 ,结合 公式 (6) 计算 
得 到 直线 ， 将 肉牛 轮廓 分 为 上 下 轮廓 ， 分 析 肉 4 
头 部 特征 确定 头 部 去 除 点 ， 如 图 5 所 示 。 首 先 ， 
通过 凸 包 分 析 方法 获取 4(x，7) 网 像 边缘 轮廓 上 


让 


WDEK, WAS (a) 所 示 ; 其 次 ， 对 主 骨 架 像 
素 点 进行 抽样 ， 并 将 其 拟 合成 直线 ， 直 线 与 轮廓 
尾部 交 于 点 P， 尾 部 凸 点 距离 P 点 的 最 近 点 即 为 


EIRA Pa” WES (b) 所 示 ; 最 后 ,根据 尾 根 

点 和 公式 (6) 计算 得 到 直线 ， 将 肉牛 整体 轮廓 

分 为 上 下 轮廓 ， 如 图 5 (c) 和 图 5 (d) 所 示 。 
为 确定 肉牛 头 部 去 除 点 ， 首 先 计算 图 6 (a) 


EDE CZ ALTS (b) 尾 根 点 


(c) 上 轮廓 d) FHA 


包 分 析 方 法 获取 肉牛 图 像 上 下 轮廓 边缘 


图 5 加 
Fig. 5 Getting the upper and lower contour edges by 
convex hull analysis method 
中 各 骨架 端点 到 主 骨架 前 端点 忆 , 的 距离 ， 最 短 
距离 为 D,,,; 其 次 在 主 骨 架 上 检索 距离 P, 点 DD, 
长 度 的 像素 点 7T,， 如 图 6 (b) 所 示 ; 最 后 ， 分 
别 遍历 肉牛 上 下 轮廓 边缘 ， 检 索 提 取得 到 上 下 轮 
廊 距 离 点 7, 最 近 的 点 O, 和 0O,， 即 为 肉牛 头 部 去 

除 点 ， 如 图 6 (c) 和 图 6 (d) 所 示 。 


(a) 主 骨架 前 端点 (b) 头 部 最 小 外 接 贺 


(c) 上 轮廓 头 部 去 除 点 (d) 下 轮 廊 头 部 去 除 点 


图 6 欧 氏 距离 确定 肉牛 头 部 去 除 点 
Fig. 6 Determining the head removal point by 


Euclidean distance 


确定 肉牛 头 部 去 除 点 后 ， 连 接 肉牛 头 部 去 除 
点 ， 得 到 不 会 头 部 的 肉牛 二 值 图 像 ， 将 二 值 图 像 
与 原 深度 图 像 进行 逮 辑 “与 ”操作 ， 得 到 不 含 肉 
牛头 部 的 深度 图 像 ， 如 图 7 所 示 

从 图 7 中 看 出 ， 本 研究 方法 sce SR 
BEES PAS SAMAR, KARE IR Tey Ay 


148 


(at ERR 


(d) 抬 头 图 像 (6) 低头 图 像 
图 7 多 姿态 下 肉牛 头 部 去 除 图 像 效 果 


Fig. 7 Head removal effect of beef cattle in multiple poses 


用 于 提取 体 尺 测 点 。 
2.4 肉牛 体 尺 测量 


基于 改进 U 弦 长 曲率 的 体 尺 测 点 提取 
通过 分 析 图 像 中 肉牛 轮廓 特征 ， 发 现 轮廓 曲 
线 中 肩 部 曲率 最 大 点 和 腹部 曲率 最 大 点 是 体 尺 测 
点 。 相 比 于 K 余 弦 曲 率 、 带 曲线 平滑 的 差分 曲率 
和 工 弦 长 曲率 ，U 弦 长 曲率 有 更 强 的 抗 旋转 性 和 
抗 噪 性， 并 且 癌 弱 长 所 对 应 的 曲率 与 离散 曲线 真 
实 曲率 之 间 存 在 着 一 种 线性 关系 “ 。 

计算 离散 曲线 的 U 弱 长 曲率 时 首先 要 确定 其 
支持 邻 域 5 。 为 防止 支持 邻 域 端点 落 在 肉牛 轮 
廊 曲 线 的 毛刺 点 上 ， 影 响 体 尺 测量 精度 ， 本 研究 
提出 了 一 种 基于 改进 U 弱 长 曲率 的 肉牛 体 尺 测 点 
提取 方法 。 

记 1= {Pi(x,, y, i=1, 2, =, n) HER 
头 部 后 肉牛 躯干 上 轮廓 曲线 ，P ,为 轮廓 曲线 一 像 
素 点 ， 如 图 8 所 示 。 

点 尸 的 支持 邻 域 见 公式 (7)。 

Q(P,) = [P4 P}] (7) 


i 


2.4.1 


puf 


图 8 Pi 点 U 弦 长 曲率 的 支持 邻 域 


Fig. 8 Support neighborhood of U chord curvature at point P, 


支持 邻 域 的 约束 条 件 见 公式 (8) 
| PPh =| P,P?|=U (8) 
本 研究 通过 分 析 大 量 数 据 ， 发 现 肉牛 轮廓 存 
在 毛刺 点 ， 主 要 分 为 四 类 类 型 ， 如 图 9 所 示 ， 以 
及 分 别 将 它们 顺 时 针 方 向 旋转 90°* 、180° 和 270? 


“i 


(b) 类 型 2 (c) 类 型 3 
图 9 肉牛 图 像 毛 刺 点 类 型 


Fig. 9 Burr point types of beef cattle image 


若 前 景点 已 和 其 支持 邻 域 前 后 端点 Pi、 P? 
落 在 图 9 所 示 的 P 点 ， 则 会 影响 U 弦 长 曲率 计算 
的 准确 性 ， 降 低 体 尺 测 点 提取 精度 。 通 过 分 析 肉 
牛 轮 廓 边缘 毛刺 点 的 特点 ， 本 研究 设计 一 个 5X5 
的 卷 积 核 ， 将 卷 积 核 分 为 正 负 区 域 ， 如 图 10 所 
示 ， 划 分 红色 区 域 为 正 区 域 ， 蓝 色 区 域 为 负 区 
域 。 利 用 卷 积 核 遍 历 肉牛 轮 廓 边缘 前 景点 和 支持 
邻 域 前 后 端点 ， 若 满足 公式 (9)， 则 认为 该 点 落 
在 毛刺 点 上 ， 删 除 该 点 ， 并 利用 公式 (10) 重新 
确定 该 前 景点 的 坐标 。 


(d) 类 型 4 


(a) 类 型 1 


enn) 


T] I | 4 
正 区 域 | | 
一 一 一 一 一 
| WN SE | ees 

P 
r] T] 
| 
[| | 
| eee | [sd 


图 10 划分 5X5 的 卷 积 核 正 负 区 域 
Fig. 10 Divide the positive and negative areas of the convo- 
lution kernel of 5X5 
N(P)=4 
S(P)= 255 
N(P) #0 
N(P) #0 
其 中 ,，N(P) 为 以 P 为 中 心 的 5Xx5 区 域内 除 
P 点 外 像素 值 为 255 点 的 个 数 ， 个 ; S(P) 为 P 点 
处 的 像素 值 ;N,(P) 为 以 P 为 中 心 5x5 的 卷 积 核 


(9) 
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正 区 域内 像素 值 为 255 点 的 个 数 ， 个 。N.(P) 为 
以 P 为 中 心 5x5 的 卷 积 核 负 区 域内 像素 值 为 255 
点 的 个 数 ， 个 。 


(10) 


由 于 肉牛 轮廓 曲线 为 离散 曲线 ， 可 能 在 肉牛 
轮廓 曲线 上 不 存在 满足 P, 点 的 支持 邻 域 约束 条 件 
的 前 后 端点 ， 因 此 利用 线性 插值 方法 来 预 估 满足 
PP 点 支持 邻 域 约束 条 件 的 前 后 端点 坐标 ， 差 值 计 
算 如 公式 (11). 

Pag ug +(l=-w PvO0su<1 (11) 

根据 | P,P! =U, RER u, Am E gR 
域 前 端点 P/， 同 理 可 得 邻 域 后 端点 Po 

U 弦 长 曲率 的 计算 公式 如 式 (12) 和 (13)。 


Ci = sji- (2) (12) 
A 


s,=sign[(x,-x!)(yf-yt)-Gf- x) 0,-91)] 
(13) 
分 别 为 点 


其 中 ， (x, yw). G4 yw 
Pi, PHI bb; D, =| PPPs |. 
本 研究 提出 的 改进 U 弦 长 曲率 算法 ， 首 先 判 


断 前 景点 是 否 为 毛刺 点 ; 若是 毛刺 点 ， 则 删除 该 
点 ， 再 利用 前 景点 周围 像素 点 钳 标 均值 重新 定位 


前 景点 的 位 置 坐标 ,来 提高 U 苞 长 曲率 计算 的 准 
确 性 。 

利用 改进 UU 弦 长 曲率 算法 提取 肉牛 轮廓 各 像 
素 点 曲率 值 ， 依 次 连接 肉牛 轮廓 各 像素 点 的 曲率 
值 ， 得 到 曲率 曲线 ， 利 用 五 点 去 差 值 法 和 线性 五 
点 平滑 法 ， 对 该 曲率 曲线 进行 平滑 处 理 ， 提 取 曲 
线 肩 部 曲率 值 极 值 点 和 腹部 曲率 值 极 值 点 ， 即 为 
体 尺 测 点 。 
2.4.2 ”肉牛 体 尺 计算 

在 传统 人 工 测量 肉牛 体 尺 时 ， 体 高 通常 指 荞 
甲 最 高 点 到 地 面 的 垂直 距离 ， 体 直 长 指 肉牛 肩 端 
点 到 坐骨 节点 的 水 平 距 离 ， 腹 宽 指 肉牛 腹部 的 最 
宽度 "9， 如 图 11 所 示 。 


图 11 肉牛 体 尺 计算 示意 图 
Fig. 11 Schematic diagram of beef cattle body 


size calculation 


Fa FC HE PY ae RE KUR he TL A E THE n “MRR 
点 ， 根 据 相 机 的 内 参 将 像素 点 转换 到 三 维 坐标 系 
下 ， 通 过 ?个 点 拟 合 出 地 平面 ， 获 取 平 面 方程 ， 
如 公式 (14); 其 次 ， 筛 选 奶牛 前 1/2 区 域内 所 有 
坐标 点 的 Z 值 数据 “"， 获 取 Z 最 小 值 点 五 (x, y) 
最 后 ， 将 最 小 值 点 转换 到 三 维 坐标 系 下 为 
H/'(x”，y，z')， 计算 五 ' 到 平面 的 距离 即 为 奶牛 体 


Fy Ho 
Ax + By +Cz+D=0 (14) 
获取 腹部 区 域 中 曲率 极 值 点 为 腹 宽 测 点 
fX Yo 2) AOs Y» z), AWA E E 


点 之 间 的 欧式 距离 即 为 腹 宽 参数 4,， 如 公 
IN CIS Je 
A,= (xi -x (y — 2h (2, - VP (15) 
当 肉牛 身体 牌 斜 严重 时 ， 肩 部 区 域 会 出 现 两 
个 曲率 极 值 点 。 为 获取 体 直 长 前 端点 ， 先 判断 肩 
部 曲率 极 值 点 个 数 ， 若 存在 两 个 极 值 点 ， 则 取 两 


个 极 值 点 的 均值 点 作为 体 直 长 前 端点 
Pi(Xies Vies Ze)» 尾 根 点 Pe Xs Yas Za AVE EL 
长 后 端点 ， 计 算 体 直 长 前 端点 到 腹 宽 中 点 
AA, mo, 15 72) 的 欧 欧式 距离 为 d,， 
gp 
计算 公式 如 (16). 
L=d,td, (16) 


3 结果 与 分 析 


头 部 去 除 结果 分 析 
随机 选取 23 头 肉 牛 的 深度 图 像 数据 ， 由 于 


150 


相 邻 帧 之 间 的 肉牛 姿态 差异 较 小 ， 利 用 抽 帧 的 方 
法 ,每 间隔 10 帧 抽取 1 帧 作为 测试 数据 ， 共 得 到 
674 帧 图 像 数 据 。 首 先 ， 通过 人 了 眼 观 察 的 方法 对 
674 帧 深度 图 像 进 行 姿 态 的 初步 分 类 ， 分 别 为 姿 
态 端 正 、 左 焉 (本 研究 中 左下 指 牛头 部 和 艇 干 不 
在 一 条 直线 上 ， 且 头颈 部 牌 向 通道 左 侧 ) 、 右 在 
(本 研究 中 右 牌 指 牛头 部 和 躯干 不 在 一 条 直线 上 ， 
且 头 颈 部 牌 向 通道 右 侧 ) 、 低 头 和 抬头 5 类 ; 然 
后 ， 再 计算 图 像 中 肉牛 姿态 稚 斜 角度 值 和 低头 抬 
头 幅 度 值 ， 来 进一步 确定 肉牛 姿态 。 

肉牛 姿态 不 端正 程度 可 以 用 0 角度 来 衡量 。 
首先 ， 提取 主 骨架 最 小 外 接 和 矩形 ， 计 算 和 矩形 Y 轴 
方向 的 中 点 所 在 X 轴 方向 直线 与 主 骨 架 交 点 M， 
确定 主 骨架 前 端点 H; 然后 ， 确 定 肩 部 对 应 骨架 
点 S; 最 后 ， 连 接 MS 和 HS, 求 直线 7, 和 17, 之 间 
的 夹 角 9， 如 图 12 所 示 。 通 过 该 方法 将 674 帧 深 
度 图 像 分 为 3 类 ， 若 |9|< 15"”"， 则 标记 该 深度 
图 为 姿态 端正 帧 ;， 若 |9|> 15 并 且 人 眼 观察 认为 
该 图 像 左 乍 ， 则 标记 该 帧 为 左 和 至 图 像 ; 4) 0|> 
15 并 且 人 眼 观察 认为 该 图 像 右 和 看， 则 标记 该 帧 
为 右 焉 图 像 。 


ji 


图 12 肉牛 姿态 不 端正 度 评价 0 角度 计算 方法 
Fig. 12 Evaluation method of beef cattle posture 
skew degree 0 

为 计算 肉牛 低头 和 抬头 的 幅度 值 ， 将 次 度 图 
像 转换 为 三 维 点 云 ， 利 用 CloudCompare 软件 标 
注 牛 只 头 部 和 肩 部 的 高 度 ， 并 计算 高 度 差 值 于 ,， 
用 五 , 值 表示 低头 和 拾 头 幅度 值 大 小 。 

以 往 研究 多 采用 凸 包 分 析 的 方法 来 去 除 图 像 
中 牛 只 头 部 "。 为 验证 本 研究 肉牛 头 部 去 除 方法 
的 高 效 性 ,分别 用 凸 包 分 析 方 法 和 本 研究 所 提 方 


法 对 图 像 中 不 同 姿态 下 的 肉牛 头 部 进行 去 除 试验 
分 析 ， 结 果 如 表 1 所 示 。 可 知 ， 本 研究 所 提 方 法 
在 不 同 姿态 下 ， 肉 牛头 部 去 除 成 功率 均 高 于 
92%， 且 均 高 于 凸 包 分 析 方 法 的 尖 部 去 除 成 
功率 。 


表 1 五 种 姿态 下 肉牛 头 部 去 除 结果 


Table 1 Head removal results of five postures 


姿态 端 LE AE 低头 Hk 
正 图 像 图 像 AR 图像 图 像 
总 帧 数 / 帧 298 194 182 252 30 
Here 126 90 36 90 18 
本 文 方法 成 功 帧 数 / 帧 282 186 168 235 30 
凸 包 分 析 方 法 成 功率 /% 42.28 46.39 19.78 35.71 60 
本 文 方法 成 功率 /% 94.63 98.88 92.31 93.25 100 


为 进一步 验证 本 研究 肉牛 头 部 去 除 方法 在 不 
同 牌 斜 角度 下 头 部 去 除 效果 ， 计 算得 376 帧 左 算 
图 像 和 右 牌 图 像 中 0 均值 为 33.67" ， 针 对 姿态 算 
斜 程度 ， 将 0 值 小 于 均值 的 图 像 标 记 为 轻 度 焉 斜 
图 像 ，0 值 大 于 等 于 均值 的 图 像 标 记 为 重度 牌 斜 
图 像 ， 并 按照 此 方法 对 左 焉 图像、 右 焉 图 像 进 行 
了 细 分 ， 细 分 后 图 像 中 肉牛 头 部 去 除 结果 如 表 2 
和 表 3 所 示 。 

表 2 左 牌 姿态 细 分 后 肉牛 头 部 去 除 结果 


Table 2 Head removal results after the left-slanted posture 


was subdivided 


AWE th LAPT ITI 本 研究 方法 成 6 均值 / 


肉牛 姿态 
数 / 帧 ”成功 帧 数 / 帧 ” 功 帧 数 / 帧 ($) 
ERER 100 62 98 21.75 
HEALER 94 28 88 47.46 


表 3 右 焉 姿态 细 分 后 肉牛 头 部 去 除 结果 


Table 3 Head removal result after the right-slanted posture 


was subdivided 


总 帧 ” 凸 包 分 析 方 法 本 研究 方法 成 


为 牛 姿态 0 均值 /(?) 

nT oum mgu E 
EREK 84 24 82 22.09 
EERE 98 10 86 42.56 


由 表 2 和 表 3 可 知 ， 肉 牛 轻 度 左 牌 图 像 9 均 
值 为 21.75°* ， 轻 度 右 焉 图 像 9 均值 为 22.09°*。 对 
于 轻 度 左 牌 图 像 和 轻 度 右 牌 图 像 中 的 头 部 ， 凸 包 
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分 析 方 法 去 除 成 功率 分 别 为 62.00% Fil 28.57%, 
本 研究 提出 的 方法 去 除 成 功率 分 别 为 98.00% 和 
97.62%。 对 于 重度 左 牌 图 像 和 重度 右 牌 图 像 中 的 
头 部 ， 凸 包 分 析 方 法 去 除 成 功率 分 别 为 29.79% 
和 10.20%， 本 研究 提出 的 方法 去 除 成 功率 分 别 
为 93.62% 和 87.76%。 上 述 分 析 表 明 ， 图 像 中 姿 
态 乍 斜 角度 90 越 大 ， 尖 部 去 除 越 困 难 ， 两 种 方法 
的 头 部 去 除 成 功率 均 降 低 ， 本 研究 方法 在 5 种 姿 
态 下 的 头 部 去 除 成 功率 均 高 于 凸 包 分 析 方法 。 
综合 考虑 图 像 中 肉牛 低头 和 抬头 幅度 以 及 姿 


80r 


低头 幅度 /cm 


ie) 
局 
T 


态 牌 斜 程度 对 本 研究 提出 的 头 部 去 除 方法 去 除 结 
果 的 有 影响， 对 244 帧 深度 图 像 的 姿态 焉 斜 角度 
909、 低头 幅度 及 ,以 及 头 部 去 除 结果 进行 统计 分 
析 ， 如 图 13 所 示 。 其 中 ， 蓝 色 点 和 绿色 点 为 肉 
牛头 部 去 除 失败 帧 ， 蓝 色 点 主要 集中 在 五, 值 较 
大 区 域 , 绿色 点 主要 集中 在 0 值 较 大 区 域 ， 因 此 
当 图 像 中 肉牛 低 尖 和 抬头 幅度 较 大 或 姿态 焉 和 斜 程 
度 较 高 时 ， 存 在 头 部 去 除 失败 的 情况 。224 WR 
度 图 像 中 仅 13 帧 头 部 去 除 失败 ， 证 明了 本 研究 
所 提 方 法 的 有 效 性 。 


一 一 低头 幅度 玉 
 。 次 态 牌 斜 程度 0 


图 13 在 不 同 低头 幅度 末 , 和 姿态 焉 斜 程度 0 下 肉牛 头 部 去 除 结果 


Fig.13 Head removal results of beef cattle under different head-down posture H, and posture skew degree 0 


3.2 体 尺 测量 结果 分 析 


为 自动 提取 肉牛 体 尺 测 点 ， 对 肉牛 轮廓 进行 
一 系列 处 理 。 首 先 ， 通 过 尾 根 点 和 主 骨 架 前 端点 
对 牛 体 进行 姿态 矫正 ;其 次 ,利用 UU 弦 长 曲率 算 
法 和 改进 品 弦 长 曲率 算法 分 别 计 算得 到 肉牛 上 下 
轮廓 的 曲率 曲线 ， 结 果 如 图 14 所 示 。 对 图 14 分 
析 可 知 ， 改 进 后 的 U 弦 长 曲率 算法 计算 得 到 的 TU 
弦 长 曲率 曲线 更 加 平滑 ， 能 更 加 准确 地 计算 出 肉 
牛 轮廓 边缘 像素 点 的 曲率 值 。 

利用 人 工 标注 方法 和 本 研究 体 尺 测 点 提取 方 
法 ,在 3.1 节 获取 测试 数据 的 基础 上 ， 针 对 每 头 
肉牛 随机 选取 2 帧 深度 图 像 ， 共 获取 到 23 头 不 同 
姿态 下 肉牛 46 帧 深度 图 像 ， 进 行 体高 、 体 直 长 


和 腹 宽 的 测量 。 选 用 点 云 处 理 软件 CloudCom- 


pare 对 俯视 点 云 进 行 体 尺 参 数 的 标注 ， 对 于 每 个 
体 尺 参数 人 工 标注 三 次 ， 求 其 平均 值 作为 真实 
值 ， 肉 牛 体 尺 参数 测量 值 与 真实 值 的 结果 见 表 4， 


肉牛 标号 方法 为 : 将 第 1 头 肉 牛 的 两 帧 图 像 分 别 
标记 为 1 号 和 2 号 , 第 2 头 肉 牛 的 两 帧 图 像 分 别 
标记 为 3 号 和 4 号 …… 第 23 头 肉牛 的 两 帧 图像 分 
别 标记 为 45 号 和 46 号 。 

由 表 4 可 知 ， 体 直 长 测量 的 平均 绝对 误差 为 
2.73 cms 体高 测量 的 平均 绝对 误差 为 2.07 om; 
腹 宽 测量 的 平均 绝对 误差 为 1.47 cm。11 号 肉牛 
体高 测量 误差 较 大 ， 误差 为 9.79 cm， 主 要 原因 
是 深度 相机 分 辨 率 较 高 ， 数 据 采 集 时 存在 少量 因 
加 工 食料 而 产生 的 粉尘 在 牛 体 上 空 漂浮 ， 体 高 测 
点 提取 方法 误 将 粉尘 到 深度 相机 平面 的 距离 看 作 
是 牛 体 背部 到 相机 平面 的 距离 ， 体 高 测 点 选取 错 
误 。44 号 肉牛 的 腹 宽 测量 误差 较 大 ， 误 差 为 
4.32 cm， 主 要 原因 是 该 深度 图 像 肉 牛 轮廓 边缘 
信息 缺失 严重 ， 本 研究 虽 重新 确定 了 为 毛刺 点 的 
前 景点 的 位 置 ， 但 仍 使 腹 宽 测 点 的 提取 出 现 较 大 
误差 ， 该 因素 也 导致 体 直 长 的 测量 误差 高 于 
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(b) 改 进 算法 后 的 吕 弦 长 曲率 曲线 


图 14 U 弦 长 曲率 算法 改进 前 后 计算 得 到 肉牛 轮 廊 曲率 曲线 


Fig. 14 Curvature curve of beef cattle contour calculated by U-chord length curvature algorithm before and after improvement 


均值 。 

为 验证 体 尺 测 点 提取 方法 针对 非 端 正 姿态 下 
肉牛 的 体 尺 测量 精度 ， 将 46 帧 深度 图 像 分 为 姿 
态 端正 图 像 和 非 端 正 姿 态 图 像 ， 分 别 计算 两 类 图 
像 的 体高 平均 误差 、 体 直 长 平均 误差 和 腹 宽 平均 
误差 ， 结 果 如 表 S$ 所 示 。 可 知 ， 本 研究 体 尺 测 点 
提取 方法 在 端正 姿态 和 非 端正 姿态 下 肉牛 体 尺 测 
量 误差 均 较 低 。 与 端正 姿态 下 的 肉牛 体 尺 测量 误 
差 相 比较 ， 非 端正 姿态 下 的 肉牛 体高 测量 误差 小 
于 0.81 cm、 体 直 长 测量 误差 高 于 0.264 cm、 腹 
宽 测量 误差 高 于 0.418 cme ZOR UIK 
曲率 的 体 尺 测 点 提取 方法 在 肉牛 的 体 尺 测量 效果 
也 较 好 。 

为 验证 体 尺 测 点 提取 方法 在 多 姿态 下 肉牛 的 
体 尺 测量 效果 ， 将 体 尺 测 点 提取 方法 测 得 的 肉牛 
体 尺 结果 与 文献 [18] 3D 特征 提取 方法 、 文 献 [19] 
基于 Mask RONN 的 体 尺 测量 方法 和 文献 [7] 基 于 
YOLOv5 的 体 尺 测 量 方法 体 尺 测量 结果 进行 对 比 
分 析 ， 且 后 者 均 是 针对 端正 姿态 下 牛 只 进行 的 体 
尺 测 量 ， 如 表 6 所 示 。 可 知 ， 本 研究 方法 测 得 肉 
牛 体 直 长 误差 与 文献 [7] 测 得 体 直 长 误差 相差 较 
小 ， 且 小 于 文献 [19] 测 得 体 直 长 误差 ， 本 研究 测 


得 肉牛 体高 误差 高 于 文献 [7] 测 得 体高 误差 
0.86%， 且 小 于 文献 [18] 和 文献 [19] 测 得 体高 误 
差 。 充 分 表明 本 研究 提出 的 体 尺 测 点 提取 方法 可 
以 在 实现 多 姿态 下 肉牛 体 尺 测量 的 同时 ， 还 能 较 
大 提高 体 尺 测量 精度 。 


4 结 论 


肉牛 体 尺 参数 可 以 作为 肉牛 生长 发 育 状况 及 
选 种 育种 重要 判断 依据 。 而 养殖 场 的 肉牛 较为 活 
跃 ， 采集 到 的 图 像 中 姿态 端正 帧 较 少 ， 量 姿态 各 
异 ， 传 统 方法 难以 测量 其 体 尺 。 本 研究 通过 分 析 
肉牛 骨架 特征 和 肉牛 轮廓 边缘 特征 ， 提 出 基于 深 
度 图 像 的 多 姿态 肉牛 体 尺 自动 测量 方法 : 

(1) 提出 了 一 种 基于 骨架 的 多 姿态 下 肉牛 头 
部 去 除 方法 。 随 机 选取 肉牛 在 不 同 姿态 下 的 深度 
图 像 674 帧 ， 并 将 其 划分 为 五 种 姿态 : 姿态 端正 
图 像 298 帧 、 左 牌 图 像 194 帧 、 右 牌 图 像 182 帧 、 
低头 图 像 252 帧 和 抬头 图 像 30 帧 。 本 研究 所 提 肉 
牛头 部 去 除 方法 在 五 种 姿态 下 的 头 部 去 除 成 功率 
分 别 为 94.63%、98.88%、92.31%、93.25% 和 
100%， 均 优 于 基于 凸 包 的 头 部 去 除 方法 。 

(2) 提出 了 一 种 基于 改进 U 弦 长 曲率 的 体 尺 
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表 4 改进 U 弦 长 曲率 算法 测量 肉牛 体 尺 参 数 结果 
Table 4 Measuring the body size parameters of beef cattle by improved U-chord length curvature algorithm 
Pn 本 文 方法 测量 值 人 工 测量 什 姿态 
体高 腹 宽 体 直 长 体高 腹 宽 体 直 长 /m 
1 1.2803 0.5991 1.2922 1.2651 0.5806 1.3146 低头 
2 1.2784 0.5968 1.2577 1.2651 0.5806 1.3146 AE 
3 1.2528 0.5661 1.2546 1.2678 0.5790 1.2615 低头 
4 1.2984 0.5968 1.2577 1.2678 0.5790 1.2615 AE 
5 1.2941 0.5852 1.2547 1.3121 0.5521 1.2520 AE 
6 1.3226 0.5761 1.2472 1.3121 0.5521 1.2520 姿态 端正 
7 1.3161 0.5531 1.2726 1.2946 0.5440 1.2260 姿态 端正 
8 1.3283 0.5595 1.2613 1.2946 0.5440 1.2260 EE 
9 1.2455 0.6524 1.3217 1.2259 0.6522 1.3523 低头 
10 1.2392 0.6706 1.3115 1.2259 0.6522 1.3523 EE 
S 11 1.1279 0.6435 1.3749 1.2258 0.6486 1.3192 姿态 端正 
©) 12 1.2420 0.6542 1.3389 1.2258 0.6486 1.3192 EE 
N 13 1.2813 0.5836 1.2308 1.2793 0.5575 1.3111 低头 
bo 14 1.3032 0.5765 1.2477 1.2893 0.5575 1.3111 RE 
© 15 1.2773 0.5829 1.1892 1.2348 0.6002 1.1899 AE 
© 16 1.2645 0.5904 1.1791 1.2348 0.6002 1.1899 姿态 端正 
N 17 1.2872 0.5736 1.1327 1.2606 0.5909 1.1365 AE 
= 18 1.2767 0.6093 1.1140 1.2606 0.5909 1.1365 低头 
N 19 1.2565 0.6164 1.1514 1.2348 0.6002 1.1899 AE 
© 20 1.2422 0.6187 1.1553 1.2348 0.6002 1.1899 HE 
N 21 1.3331 0.5295 1.2316 1.3187 0.5189 1.2039 姿态 端正 
ni 22 1.3414 0.5291 1.2271 1.3187 0.5189 1.2039 AE 
2 23 1.2987 0.5972 1.2444 1.2554 0.6207 1.2914 AE 
X< 24 1.2768 0.6066 1.2731 1.2554 0.6207 1.2914 低头 
© 25 1.2527 0.6408 1.2174 1.2635 0.6209 1.2267 低头 
= 26 1.2793 0.6392 1.1921 1.2635 0.6209 1.2267 AE 
T 27 1.2690 0.6048 1.2176 1.2506 0.6144 1.2206 E 
O 28 1.2631 0.5996 1.2079 1.2506 0.6144 1.2206 姿态 端正 
29 1.3248 0.5848 1.2481 1.3444 0.5956 1.2821 端正 
30 1.3271 0.6140 1.3003 1.3444 0.5956 1.2821 低头 
31 1.2858 0.6969 1.1688 1.2549 0.6986 1.2005 AE 
32 1.2360 0.6890 1.1303 1.2549 0.6986 1.2005 E 
33 1.2803 0.6457 1.2163 1.2631 0.6465 1.2443 姿态 端正 
34 1.2509 0.6486 1.2287 1.2631 0.6465 1.2443 EE 
35 1.3167 0.5935 1.3026 1.3263 0.5825 1.3017 低头 
36 1.3260 0.5947 1.3298 1.3263 0.5825 1.3017 端正 
37 1.2250 0.6359 1.2869 1.2403 0.6508 1.2896 低头 
38 1.1959 0.6546 1.2919 1.2403 0.6508 1.2896 端正 
39 1.2190 0.5985 1.2307 1.2321 0.6069 1.2356 AE 
40 1.1941 0.5960 1.2789 1.2321 0.6069 1.2356 端正 
41 1.3209 0.6147 1.1853 1.3187 0.6285 1.2179 AE 
42 1.3093 0.6183 1.1557 1.3187 0.6285 1.2179 AE 
43 1.2517 0.6017 1.1886 1.2674 0.6176 1.2035 姿态 端正 
44 1.2472 0.5744 1.1662 1.2674 0.6176 1.2035 AE 
45 1.4028 0.5970 1.2225 1.3737 0.6276 1.2762 AE 
46 1.3898 0.6091 1.2561 1.3737 0.6276 1.2762 AE 
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表 5 不 同 姿态 下 肉牛 体 尺 测量 平均 误差 
Table 5 Average errors of body measure measurement under 


different postures 


肉牛 姿态 体高 平均 体 直 长 平均 宽 平均 
误差 /cm 误差 /cm 误差 /cm 
端正 姿态 2.742 2.515 1.126 
非 端正 姿态 1.932 2.779 1.544 


表 6 不 同方 法 肉牛 体 尺 测量 平均 相对 误差 


Table 6 Average relative errors of body measurement under 


different methods 
方法 WAWE 体 直 长 所 
误差 /% 
3D 特征 提取 方法 <3.00 = 
REF Mask RCNN HERMETE 5.78 6.09 
基于 YOLOv5 的 体 尺 测量 方法 "| 0.76 2.14 
改进 U 弱 长 曲率 方法 1.62 2.18 


测 点 提取 方法 。 随 机 选取 23 头 肉 牛 在 不 同 姿态 
下 的 深度 图 像 共 46 帧 ， 利 用 改进 U 弱 长 曲率 算 
法 提取 体 尺 测 点 ， 并 计算 体 尺 参 数 ， 测 得 肉牛 体 
直 长 平均 绝对 误差 为 2.73 cm， 体 高 平均 绝对 误 
差 为 2.07 cm， 腹 宽 平均 绝对 误差 为 1.47 cm， 且 
本 研究 所 提 方 法 测 得 非 端 正 姿态 下 肉牛 图 像 的 体 
直 长 平均 绝对 误差 为 2.779 cm， 体 高 平均 绝对 误 
差 为 1.932 cm， 腹 宽 平 均 绝 对 误差 为 1.544 cmo 
充分 表明 ,改进 U 弱 长 曲率 的 体 尺 测 点 提取 方法 
可 以 提取 不 同 姿态 下 肉牛 体 尺 测 点 ， 且 在 不 同 姿 
态 下 体 尺 测量 误差 均 较 小 。 
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Automatic Measurement of Multi-Posture Beef Cattle 
Body Size Based on Depth Image 


YE Wenshuai”, KANG Xi”, HE Zhijiang™, LI Mengfei'”, LIU Gang™ 


(1. Key Lab of Smart Agriculture Systems, Ministry of Education, China Agricultural University, Beijing 100083, 
China; 2. Key Laboratory of Agricultural Information Acquisition Technology, Ministry of Agriculture and Rural Af- 
fairs,China Agricultural University, Beijing 100083, China; 3. School of Computing and Data Engineering, Ning- 
boTech University, Ningbo 315200, China) 


Abstract: Beef cattle in the farm are active, which leads the collection of posture of the beef cattle changeable, so it is difficult to 
automatically measure the body size of the beef cattle. Aiming at the above problems, an automatic measurement method for 
beef cattle's body size under multi-pose was proposed by analyzing the skeleton features of beef cattle head and the edge con- 
tour features of beef cattle images. Firstly, the consumer-grade depth camera Azure Kinect DK was used to collect the top-view 
depth video data directly above the beef cattle and the video data were divided into frames to obtain the original depth image. 
Secondly, the original depth image was processed by shadow interpolation, normalization, image segmentation and connected 
domain to remove the complex background and obtain the target image containing only beef cattle. Thirdly, the Zhang-Suen al- 
gorithm was used to extract the beef cattle skeleton of the target image, and calculated the intersection points and endpoints of 
the skeleton, so as to analyze the characteristics of the beef cattle head to determine the head removal point , and to remove the 
beef cattle head information from the image. Finally, the curvature curve of the beef cattle profile was obtained by the improved 
U-chord curvature method. The body measurement points were determined according to the curvature value and converted into 
three-dimensional spaces to calculate the body size parameters. In this paper, the postures of beef cattle, which were analyzed 
by a large amount of depth image data, were divided into left crooked, right crooked, correct posture, head down and head up, 
respectively. The test results showed that the head removal method proposed based on the skeleton in multiple postures hads 
head removel success rate higher than 92% in the five postures. Using the body measurement point extraction method based on 
the improved U-chord curvature proposed, the average absolute error of body length measurement was 2.73 cm, the average ab- 
solute error of body height measurement was 2.07 cm, and the average absolute error of belly width measurement was 1.47 cm. 
The method provides a better way to achieve the automatic measurement of beef cattle body size in multiple poses. 
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